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I. INTRODUCTION 

Control and optimization of quantum systems have 
been recognized as important issues for many years pi 
and control theory for quantum systems has been de- 
veloped since the 1980s [2H1]- There has been consider- 
able recent progress in both theory and experiment E] ■ 
However, despite this progress, there are still many chal- 
lenges. Most quantum control schemes rely on open-loop 
control design based on mathematical models of the sys- 
tem to be controlled. However, accurate models are not 
often not available, especially for manufactured quan- 
tum systems such as artificial quantum dot atoms or 
molecules. Therefore, system identification [7| is a crucial 
prerequisite for quantum control. 

In the quantum information domain, procedures for 
characterization of quantum dynamical maps are often 
known as quantum-process tomography (QPT) [SlllOj 
and many schemes have been proposed to identify the 
unitary (or completely positive) processes, for exam- 
ple, standard quantum-process tomography (SQPT) [TT]- 
[13], ancilla-assisted process tomography (AAPT) [l4l - 
116] and direct characterization of quantum dynamics 
(DCQD) [T7]. However, if control of the system's dynam- 
ics is the objective, what we really need to characterize is 
not a global process but the generators of the dynamical 
evolution such as the Hamiltonian and dissipation oper- 
ators. The problem of Hamiltonian tomography (HT), 
though less well-understood, has also begun to be con- 
sidered recently by a few authors [T8H2T] . Although QPT 
and HT differ in various regards, both try to infer infor- 
mation about the quantum dynamics from experiments 
performed on systems, and both can be studied from the 
point of view of system identification with broad tasks in- 
cluding (1) experimental design and data gathering, (2) 
choice of model sets and model calculation, and (3) model 
validation. 

Recently the quantum system identification problem 
has been briefiy explored from cybernetical point of view, 
and underlining the important role of experimental de- 
sign [22] • In this article we follow this line of inquiry. 
Throughout the paper, we make the following basic as- 



sumptions: (1) the quantum system can be repeatedly 
initialized in a (limited) set of known states; (2) that we 
can let the system evolve for a desired time t; and (3) that 
some projective measurements can be performed on the 
quantum system. The main question we are interested 
in in this context is how the choice of the initialization 
and measurement affect the amount of information we 
can acquire about the dynamics of the system. Given 
any a limited range of options for the experimental de- 
sign, e.g., a range of measurements we could perform, 
different choices for the initial states, or different control 
Hamiltonians, how to choose the best experimental de- 
sign, and what are the theoretical limitations? Finally, 
we are interested in efficient ways to extracting the rele- 
vant information from noisy experimental data. 

The paper is organized as follows: In Sec. II we dis- 
cuss the model and basic design assumptions. Sec HI 
deals with the general question of model identifiability in 
various settings, and in Sec IV we compare several differ- 
ent stategies for parameter estimation from a limited set 
of noisy data from simulated experiments see how they 
measure up. 

II. MODEL AND DESIGN ASSUMPTIONS 

To keep the analysis tractable we consider a simple 
model of a qubit subject to a Hamiltonian H and a 
system-bath interaction modelled by a single Lindblad 
operator V , i.e., with system dynamics governed by the 
master equation 

M^) =_|[7j,p]+i?[y](p), (1) 

where the Lindbladian dissipation term is given by 

V[V]{p)^VpV^ - \{VW + VV^). (2) 

We shall further simplify the problem by assuming that 
is a Hermitian operator representing a depolarizing 
channel or pure phase relaxation in some basis. Without 
loss of generality we can choose the basis so that V is 
diagonal, in fact we can choose V = \/^crz with ~ 
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diag(l,— 1) and 7 > 0. Under these assumptions the 
master equation simphfies 

V[a,]{p) = l{a,pa,-p). (3) 

The control Hamiltonian can be expanded with respect 
to the PauU basis {cr^, CTj,, cr^} 

H{t) = '^(iOQ{t)a,+^i{t)a.,~uj2{t)uy) (4) 

with possibly time-dependent coefBcients Wq, {t) . It is con- 
venient to consider a real representation of the system. 
Following the approach in |23| we expand p with respect 
to the standard Pauli basis for the 2x2 Hermitian ma- 
trices 

P^\[I + Vx(Jx+Vy(7y+V^(7z), (5) 

where the coefhcients are Va = Tr(/5CTQ,). Similarly ex- 
panding the dynamical operators allows us to recast 
Eq. ([I]) in following Bloch equation [h = 1) 

Vx{t)\ ( -7 -^o(i) -W2(t)\ (vx{t)\ 
Vy{t) = Mt) -1 -Mt) {vvit)]- (6) 

Vzit)J \iO2it) OJl{t) / \v,{t)J 

Using this simple model for illustration we subse- 
quently consider the experimental design from three as- 
pects: (1) initialization procedures, (2) measurement 
choice and (3) Hamiltonian design. 

(1) Initialization. We assume the ability to prepare 
the system in some initial state 

IV-zW) =cosf |0)-Hsinf |1), (7) 

with respect to the basis {|0), |1)}, which coincides with 
the eigenbasis of V. We can formally represent the ini- 
tialization procedure by the operator 11(07), which is the 
projector onto the state with / indicating initializa- 
tion. With these restrictions the design of the initializa- 
tion procedure is reduced to the selection of parameter 
9i. Note that we assume that we can only prepare one 
fixed initial state, not a full set of basis states. 

(2) Measurement. We assume the ability to perform 
a two-outcome projective measurement 

M = M+ - Af_ |m+> (m+ 1 - |m_) (m„ | , (8) 

where the measurement basis states can be written as 

|to+) = cos^lO) -f sin%L|i) (9a) 
|m_) =sin^|0) -cos^ll) (9b) 

so that the choice of the measurement can be reduced to 
suitable choice of the parameter 9m , and we shall indicate 
this by writing M{9m)- 

(3) Hamiltonian. In practice we may or may not 
have the freedom to choose the type of Hamiltonian but it 
will be instructive to consider the identification problem 
for the following three cases: 



(a) uJz{t) = > and LiJx{t) = ujy{t) = 0. 

(b) ijJx{t) = > and ijJy{t) — ujz{t) = 0. 

(c) LiJy{t) = cjj, > and w^(t) = w^(<) = 0. 

In this context the experimental design problem can 
be reduced to the problem of choosing suitable values 9i 
and 9m to identify the system parameters 7 and in 
each of the three cases above. The problem considered 
here is similar to that considered in |21j . in particular 
we still only allow a single initial state and single fixed 
measurement, but unlike in |21| the initial state and the 
measurement are not assumed to commute with the de- 
phasing operator. 

III. MODEL IDENTIFIABILITY 

We now consider the identification problem for the de- 
phasing qubit system, concentrating on experimental de- 
sign issues for simultaneously discriminating dephasing 
parameter and Hamiltonian parameter of two-level de- 
phasing systems. 

A. H = uj^a^, V = V^CTz 

In this special case H commutes with the dephasing 
operator V. Solving the equation 

Vx{t)\ (-1 -^z 0\ (Vx{t)\ 

iy{t) U -7 Vy{t) (10) 

Vz{t)J \ 0/ \vz{t)J 

with the initial state ([t]) gives 

Vx{t)\ /e"'''* cosWj;i sinfl/X 
Vy{t) = e-^^sinujztsin9j (11) 
Vz{t)J V cos 9 J J 

and applying the binary-outcome projective measure- 
ment M{9m) yields the measurement traces p±{t) = 
Tr[A/±p(t)]. Defining p±{t) = 2p±{t) - 1, we have 

p±{t) = ± COS0/ cos0M±e~'''* cosWzisinS/ sin 6*71/, (12) 

from which we can see that we can obtain information 
about the system parameters ujz and 7 if and only if 
sm9j 7^ and sin^^f 7^ 0, i.e., if neither the initial 
state preparation H(0/) nor the measurement M com- 
mutes with H and V. The measurement traces also yield 
information about 9i and 9m, i-e., we can determine the 
relative angles between the initialization and measure- 
ment axis and the fixed Hamiltonian/dephasing axis, if 
they are not known a-priori. We also see that the vis- 
ibility is maximized if sin 0/ sin — 1, which will be 
the case if the initialization and measurement axis are 
orthogonal to the joint Hamiltonian and dephasing axis. 
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These results make physical sense. As [H, V] = 0, if 
the initial state preparation li{0i) commutes with H and 

V then the initial state is a stationary state of the dy- 
namics and the measurement outcome is constant in time 
c± = 5(1 ± !)• If sin 6*/ ^ then the initial state is not 
stationary and the state follows a spiral path towards the 
joint Hamiltonian and dephasing axis but as both H and 

V are proportional to Cz, Tr[(Tzp(t)] is a conserved quan- 
tity of the dynamics. If sin 9m = then the measurement 
commutes with cr^, and as Tr[(Tzp(t)] is a conserved quan- 
tity, we again obtain no information about the dynamics 
from which to identify the system parameters. 



B. H = ujxOx, V 



In this case, the equation (|6| is reduced into the fol- 
lowing equation 



(13) 





and the solution of Eq. ( 13 1 for the initial state Q is 

e^'*'* sin Qi \ 

$f(t)cos0/ (14) 
,$|(t)cos0j 



where we set = 




$^(t) = -e~2*g^sinw,i 
7 

$i^(t) = 2*[cosw^t -I- sinwj.^]. 



(15a) 
(15b) 



If < 7^/4 then Q will be purely imaginary and the 
sine and cosine terms above are replaced by the respec- 
tive hyperbolic functions. If = 7^/4, the expression 
sin(u;t) must be analytically continued. 
Eq. ( 14 ) implies that we can obtain the estimated value 



of the following probabilities 
p±{t) = ±6^''* sinOi sm0M ± <^3{t) cosOi cosOm- (16) 

Noting that $f (t) depends on both ojx and 7 this shows 
that we can obtain full information about the system 
parameters if and only if cos 6*/ ^ and cos 9m 1^ 0. If 
cos 9i — or cos 9m = then we can still identify 7 but 
not ujx- Again this makes sense. 

cos 9i = for 6*/ = |, i.e., if the initial state is an 
eigenstate of the Hamiltonian. Since [H, V] in this 
case, eigenstates of H are not stationary. However, since 
the Hamiltonian and dephasing axis are orthogonal, the 
initial state remains in a plane orthogonal to the dephas- 
ing axis, the z — plane in our case, following the path 
x{t) = e^'*'*. Thus we have [H,p{t)] = for all times, and 
we can therefore not obtain any information about the 
Hamiltonian parameter lUx , but we can still obtain infor- 
mation about the dephasing parameter 7. If the Hamil- 
tonian and dephasing axis were not orthogonal then we 



would expect to be able to identify both the Hamiltonian 
and dephasing parameters even if the initial state was an 
eigenstate of H as in this case it would not remain an 
eigenstate of H under the evolution. 

If cos 9 1 ^ but cos 9m — then the measurement 
commutes with the Hamiltonian. Transforming to the 
Heisenberg picture, 

M±{t) = ~i[M±{t),H] - ^V[a,]M±{t), 

and again one can show that M±{t) remains orthogonal 
to the dephasing axis and Tr[M±(i)/9o] = Tr[M±p(i)] is 
independent of the Hamiltonian H, explaining why we 
cannot obtain any information about H in this case. 



In this case equation ^ is reduced into the following 
equation 





(17) 



whose solution for the initial state Q is 




f^lsit) sin 61/ - 
\n,{t)cos9j 



2 ^ sint2ii,tcos6'/ 




e~2' 



smCjytsm9i , 



where ujy = 



and 



e 2*[cosajj^t- 



sincl;,,tl 



$^3 {t) — e 2 * [cos ujyt + 2^ sin coyt] 



(18) 

(19a) 
(19b) 



This implies that we can obtain the estimated value of 
the probabilities 

7 7 
p±{t) = ±aie~2* cosQyt ± Q;2e" 2* sinujyt (20) 

where the coefficient functions are 



ai 
a2 



COS(f/ 

7 



M) 

cos(9i + 9m) 



7M 



(21a) 
(21b) 



As before, if < 7^/4 then w will be purely imagi- 
nary and the sine and cosine terms above turn into their 
respective hyperbolic sine and cosine equivalents, and if 
ujy = 7^/4, the expression lj^^ s\n(uit) must be analyti- 
cally continued. 

In this case it is quite interesting to notice that it is 
impossible to find such 9i and 9m that cos (6*/ — 9m) = 
su\{9i — 9m) = cos{9i + 9m) — and thus that we can 
identify both model parameters for any choice of the ini- 
tial state and measurement. This also makes sense be- 
cause regardless of the choice of 9i and 9m the initial 
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FIG. 1: Evolution of system state on the Bloch sphere and projection onto measurement axis for H — V: In the pathological 
case when the measurement axis coincides with the H, l^-axis, no information about the system parameters can be obtained 
(left); Otherwise the measurement trace contains information about both H and V (right), provided the initial state is not 
stationary. 
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FIG. 2: Evolution of system state on the Bloch sphere and projection onto measurement axis for // _L V: If the measurement 
axis coincides with the H-axis, only information about the decoherence parameter 7 can be obtained (left). In case C the 
measurement trace always contains information about both H and V (right) for any initial state and measurement. 



state in this case is always orthgonal to the Hamiltonian 
axis, and the measurement M is always orthogonal to 
H, [M,H] ^ 0, there are no conserved quantities and 
the only stationary state of the system is the completely 
mixed state. 



IV. PARAMETER ESTIMATION 

In this section we explore how to estimate dephasing 
and Hamiltonian parameters from limited noisy measure- 
ment data. In an actual experiment we can only estimate 
the probabilities p± (t) at a finite number of times tk by 



repeatedly initializing the system in some fixed state po, 
and letting it evolve for time tk before performing the 
projective measurement M. Each single repetition of the 
experiment yields a binary outcome +1 or —1 and we 
can estimate the probability p±{tk) by repeating the ex- 
periment TVe times and computing the relative frequen- 
cies of the respective measurement outcomes ±1, e.g., 
P±itk) = ^■ 

To model noisy experimental data we could generate 
'P±it) by adding a zero-mean White- Gaussian noise sig- 
nal gt to p±{t), i.e., p±{t) = p±it) + gt- By the Law 
of Large Numbers and Iterated-logarithm Law [23] this 
gives a Gaussian distribution p± (t) with mean p± (t) and 
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variance ^ ^°^2N ^'^^ ^ oo. For large this 
should be a good error model, but for small it may 
not accurately capture the nature of the projection noise, 
which follows a Poisson distribution. To more accurately 
model noisy experimental data when the number of mea- 
surement repetitions is relatively small we can simulate 
the actual experiment by generating random numbers 
r„ between and 1, drawn from a uniform distribution, 
and setting p+ — No/Ne, where A'o is the number of 
rn < P+- 



A. Fourier and time-series analysis 

One commonly used technique to find frequency com- 
ponents in a noisy time-domain signal is the Fourier 
transform or its discrete version, the discrete or Fast 
Fourier Transform (FFT). In principle, this allows us to 
estimate the Hamiltonian parameters (frequencies) w^,. 



and I 



and the damping rate 7 can be estimated 



from the Lorentzian broadening of the Fourier peak |21j . 
When the signal is sparse, i.e., we only have a relatively 
small number of sample points, and noisy, however, this 
approach becomes problematic. Fig.[3]shows that we can 
still identify a peak in the spectrum for sparse noisy sig- 
nals but the accuracy of the estimate is limited, and accu- 
rately estimating the dephasing rate 7 from the broaden- 
ing of the peak, given the distortion, is very challenging. 

Alternatively, once an estimated value for has been 
obtained, 7 can be deduced in other ways. In the simplest 
case A, where the evolution is given by Eq. ( 12 ), we have 



cos ujzt sin 6 J s'm 6m ^ P+{t) — cos 6*/ cos 6*^/ (22) 



which we can rewrite as 

^ _ jQg f P+(t)-COS 81 COS Om 

' ^ \ COS cj^t sin 1 sin 6 M 



A 



resolution. The results can be noticeably improved by 
averaging over multiple time-series of p+(i). If the mean 
value of the 7-estimates over several time series is used 
as the final estimate for 7, the estimated error of 7 can 
be remarkably reduced as shown in Fig. [s] (right), which 
indicates the mean value of 7 with 1 — 40 time-series re- 
spectively. 

The number of simulated time-series of p+ (t) required 
is determined by the desired estimation accuracy. To 
achieve a mean value of the distribution of ca. 50 with 
standard deviation less than 1.0 simulations suggest that 
at least 5 time-series with iVe = 100 are necessary, i.e., 
about 500 measurements have to be performed at each 
sample time t to ensure that the standard deviation of 
the estimated error is less than 1.0. Instead of five time 
series with N^, = 100 we could choose four times series 
with Ne — 125 or two with — 250 to reach the nearly 
same estimation accuracy, as shown in Fig. |6] In the ac- 
tual experiment one can obtain the mean-value sequence 
of 7 by continually updating p{t) at each sampled time 
t. If the mean-value sequence is found to converge to 
a certain fixed value and its standard deviation satisfies 
the requirement, the experiment measurement should be 
stopped. 

To estimate 7 by finite measurement data one will have 
to identify a suitable time when to terminate the ex- 
periment, ideally when the true value of e"'''*'* is almost 
0. A p+(t) is almost equal to P-(t) for t > t^, the sig- 
nal after this time will be effectively a pure noise signal. 
Therefore increasing the signal length beyond a certain 
critical time tkg will not improve the estimation accuracy. 

Furthermore, we consider the effect of different uJz and 
7 on the estimation accuracy. If the true value of 7 is less 
than the value of 50 above, the estimation of ujz and 7 will 
be more accurate (as shown in Fig. [?]) as the exponential 
z{t). (23) e"'*'* decays more slowly. If the value of Wz is reduced. 



Given p+{t), if Oj, 6m and uJz are known, we can there- 
fore in principle calculate z{t) and thus 7. However, in 
practice the problem is that p+{t) and all the other pa- 
rameters are not known precisely. In this case we could 
estimate 7 using the Least Square Method [25] and the 
principle of Sequential Analysis [17] . 

To explore this approach we generated simulated noisy 
data signals p+{t) with standard deviation Cf accord- 
ing to the procedure above, assuming 61 — 6m — f , 
UJz — 50, 7 = 50 and the sampled period — 1.4/is. 
If the number of measurements Ng at each time t is 
Ne = 100 then the standard deviation of estimated error 

The original noiseless probabil- 



is 



log log Nc 

2NZ ~ 



0.0874. 



ity p_|-(t), the noisy signal with standard deviation 0.0874 
and the power spectral density are shown in Fig(4j Again 
the frequency can be easily determined by the peak 
value in Fig. [ijjright). We attempt to determine 7 based 
on p+(t) using simulated time-series. For a single time 
series p+ (t) the distribution of 7 around the true value of 
50 is shown in Fig. [5]jleft). The estimation is unsatisfac- 
tory even if the signal is sampled densely at a high time 



i.e., the oscillation period is increased while its envelope 
remains unchanged the accuracy of the estimation is not 
influenced (as shown in Fig. [s]). 



B. Bayesian Analysis of Sparse Signals 

The previous section shows that we can in principle 
simultaneously identify both dephasing parameter and 
Hamiltonian parameters for simple open systems using 
Fourier and time-series analysis. Both Fourier and time- 
series analysis can deal with very noisy signals but they 
require rather dense sampling of the signal with is both 
time-consuming and resource-intensive. Ideally we would 
like to be able to estimate the parameters by sparse sam- 
pling the signal. One promising approach in this regard is 
Bayesian estimation. The main idea behind the Bayesian 
approach is to choose the parameters to be determined, 
here w and 7, to maximize a certain likelihood function 



L(p,d,CT) 



a exp 



|P d 



2cr2 



(24) 



6 
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time (arb. units) m (arb. units) 



FIG. 3: Ideal signal and sparse noisy data points for a simulated experiment assuming Model A with lo = 1 and 7 = 0.1 (left). 
We can still identify a peak in the Fourier spectrum (right) around lu = 1 but the accuracy of this estimate is limited, and 
estimating 7 from the peak broadening would be a challenge. 
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FIG. 4: Original time domain signal and simulated noisy signal (left) and power spectral density of noisy signal (right). 



where d is the measured data and the p are the prob- 
abilities predicted by the model, which depend on the 
parameters to be determined, and a is the error variance. 

Following the same approach as in we write the 
signals as a linear combination a small number of basis 
functions determined by the functional form of the sig- 
nals. Here the measurement signals p± (t) can be written 
as a linear combination of two basis functions 

P±it) = aigiit)+a2g2it). (25) 

The values of the basis functions and the coefficients are 
given in Table |T] 

Figs [9HT1I show that Bayesian analysis allows us to es- 
timate the model parameters not only for Case A but for 





Model A 


Model B 


Model C 


51 (i) 

52 (t) 
ai 

0-2 


1 

e~^' cos(a;zi) 
cos{9i) cos{9m) 
sm{9i) sm{9M) 


sm(9i) sin{9M) 
cos{9i) cos{9m) 


g-7t/ 

( 
( 


^c 
^s 
21 
IT 


os{ujt) 
in(cjf) 
i) 
^) 



TABLE I: Chosen basis functions and coefficients for Bayesian 
estimation. s{t) = cos{Qt) + ^ sin(tDt). uj — — (f )'^. 



all cases, even if the signal is very noisy, the sampling 
is sparse and the measurement and initialization choices 
are not ideal to maximize the visibility or signal-to-noise 
ratio such as Bm = j and 6j — ^. The squeezed peaks 
indicate that the accuracy of the frequency estimation is 
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Mean value of distribution is 49.7423 
Standard deviation of distribution Is 2.11 04 
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FIG. 5: Distribution of 7 based on one simulated time series (left) and mean value of 7 for 1 — 40 simulated time series (right). 



Mean value of distribution Is 49.9761 
Standard deviation of distribution Is 0.8757 
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Mean value of distribution is 49.91 1 1 
Standard deviation of distribution is 0.8854 
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FIG. 6: Distribution of mean value of 7 for Ne = 125 using 4 simulated time series (left), and for Ne = 250 using 2 simulated 
time series (right). 





Ql 

act est unccrt. 


ai 

act est unccrt. 


Model A 
Model B 
Model C 


0.3536 0.3536 0.0133 
0.6124 0.5952 0.0491 
0.9659 1.0066 0.0594 


0.6124 0.5267 0.1393 
0.3536 0.3531 0.0514 
0.2332 0.2506 0.0573 



TABLE II: Actual and estimated values of the linear coeffi- 
cients ai and Q2. 



much higher than the accuracy of the 7-estimates. 

An additional advantage of the Baycsian estimation 
does not require a-priori knowledge of the initialization or 
measurement angles Oj and 6m- Rather, the estimation 
procedure provides values for the coefficients of the basis 
functions, which are related to the parameters 61 and 0m- 



V. CONCLUDING DISCUSSION 

We have studied the issue of model identifiability and 

experiment design for open system dynamics for a de- 
phasing qubit. From the examples in Sec. Ill we can 
derive some general insights about the limits of identi- 
fiability and the role of experiment design. Unlike for 
process tomography, where we require the ability to pre- 
pare the system in many different input states, and the 
ability to measure a complete set of observables, we can 
in general extract information about both the Hamilto- 
nian and the dephasing parameters by repeating a single 
experiment, initializing the system in single fixed initial 
state and measuring a fixed basis. There are certain limi- 
tations, however. We gain no information about the value 
of the system parameters if the initial state is a stationary 
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Appendix A: Estimating peak frequency and 
dephasing rate from Fourier spectrum 

Consider a measurement trace of the form 



a + be sin(a;oi), 



(Al) 



which corresponds directly to the expected result for 

Model A if we set a = cos (6*/) cos (6'm) and b = 
sin(6'/) sin(0M)- If 0i = &m = tt/S then a = and 6 = 1. 
A similar analysis can be done for Models B and C. 

The Fourier transform of u{t)e~^* sm^ujot), where u{t) 
is the Heavyside function, is given by F{ui) 
We are interested in its absolute value 



state of the system, or if the measurement is a conserved 
quantity, although knowledge of the stationary states or 
conserved quantities restricts the dynamics and thus pro- 
vides indirect information about the system. Even if the 
initial state is not stationary or the measurement is not 
a conserved quantity, we may fail to obtain information 
about the Hamiltonian parameters is if the operators H 
and V arc orthogonal and M commutes with H, for in- 
stance. These limitations also apply to higher dimen- 
sional systems, although for such systems additional re- 
strictions on the identifiability of model parameters may 
arise. 

If the experiment design is such that the model param- 
eters are identifiable there axe various ways to extract the 
relevant parameters from a set of noisy samples (time se- 
ries), including Fourier analysis, time series analysis and 
Bayesian estimation. Although all of these approaches 
are in principle able to provide the required information, 
Bayesian estimation appears to be superior to the alter- 
natives, in particular when we are dealing with a limited 
number of noisy data points. One reason the Bayesian 
estimation is capable of providing far more accurate es- 
timates for the parameters given the same input data is 
that it utilizes information about the structure of the sig- 
nal in the form of the choice of the basis functions we are 
projecting onto. This allows us to overcome restrictions 
on the uncertainties of the parameter estimates imposed 
by Nyquist's law. 



Wo 



[72 + (a;o-w)'][7' + (w + a;o)'] 



(A2) 



Differentiating with respect to w and setting the numer- 
ator to shows that has extrema for u){u>q — 
7^ — w^) — 0, and in particular we have a maximum 
at w = — 7^ with peak value (27)"-^. 

Thus, we could in principle estimate both the fre- 
quency Wo and dephasing rate 7 from the peak height 
and position w*, 7 = (2|F|*)~-^ and ojq = sjw\ -|- 7^, but 
in practice estimating the height of the peak is delicate 
and this approach is usually very inaccuate. 

We can get a better estimate for 7 using the width of 
the peak. Let a;i,2 be the (positive) frequencies for which 
[^^(a;)! assumes half its maximum or 1/(47). Then the 
full- width-half-maximum 2d of |-F(w)| is |w2 — w\ \ or 



- 72 + 2\/3( 



^0^2 + 2^37^^2+72-^, 



Given the location and half-width d of the peak we 
can solve this equation for 7 



7 = ^-v^65(a;,,d) - 18^2 



(A3) 
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where g(w^,d) = y'Qw^ + 12^2^2 _^ Yl(p>^^^ + 3^4. Thus, 
in principle we can determine both the frequency and the 
dephasing rate by estimating the position and width of 
the fourier peak. 
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FIG. 8: Distribution of mean value of 7 for = 50Hz, 7 = 20 (left) and uj^ = 20Hz, 7 = 50 (right). 




[1] A. Blaquiere, S. Diner, and G. Lochak, (edit) Information 
Complexity and Control in Quantum Physics, Springer- 
Verlag, New York, 1987. 

[2] G. M. Huang, T. J. Tarn, and J. W. Clark, On the Con- 
trollability of Quantum Mechanical Systems, J. Math. 
Phys., vol. 24, pp. 2608-2618, 1983. 

[3] C. K. Ong, G. Huang, T. J. Tarn, and J. W. Clark, 
Invertibility of Quantum Mechanical Control Systems, 
Math. Sys. Theor., vol.17, pp. 335-350, 1984. 

[4] J. W. Clark, C. Ong, T. J. Tarn, and G. M. Huang, Quan- 
tum non-Demolition Filters, Math. Sys. Theor., vol.18, 
pp.33-53, 1985. 

[5] D. D'Alessandro, Introduction to quantum control and 
dynamics, CRC Press, 2007. 

[6] H. Rabitz, R. de Vivie-Riedle, M. Motzkus and K. 
Kompa, Whither the Future of Controlling Quantum 



Phenomena? Science vol.288, 824, 2000. 
[7] L. Ljung, System Identification: Theory for the User, 
Prentice Hall, Upper Saddle River, New Jersey (2nd Edi- 
tion) 1999. 

[8] G. M. D'Ariano, M. G. A. Paris, and M. F. Sacchi, Adv. 

Imaging Electron Phys. 128, 205 (2003). 
[9] G. M. D'Ariano and P. Lo Presti, in Quantum State Es- 
timation, edited by M. Paris and J. Rehacek, Lecture 
Notes in Physics Vol. 649 Springer, Berlin, 2004 , p. 297. 

[10] L. M. Artiles, R. D. Gill, and M. I. Gu, An invitation 
to quantum tomography, J. R. Stat. Soc. Ser. B (Stat. 
Methodol.) vol.67, 109 (2005). 

[11] M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information, Cambridge University Press, 
Cambridge, UK, (2000) 

[12] I. L. Chuang and M. A. Nielsen, Prescription for ex- 



10 



1.2 
1 

M 0.8 
c 

O) 

'cn 

^ 0.6 

Q) 

E 

§ 0.4 

0) 
0) 

E 0.2 


-0.2 




10 15 20 

time (arb. units) 



25 



30 




0.05 0.1 



0.2 0.25 0.3 0.35 0.4 

y 



FIG. 10: Ideal signal and sparse noisy sampling {N = 100, T « 25, A^e ~ 100) (left) and corresponding log-likelihood function 
(right) for model system B with uj^ = i and 7 = 0.1. 




10 15 20 
time (arb. units) 




0.05 0.1 



0.2 0.25 0.3 0.35 0.4 

y 



FIG. 11: Ideal signal and sparse noisy sampling (A'' = 100, T 
model system C with u)y = 1 and 7 = 0.1. 



25, Ne = 100) (left) and log-likehhood function (right) for 



perimental determination of the dynamics of a quantum 

black box, J. Mod. Opt. vol.44, 2455 (1997). [18] 

[13] J. F. Poyatos, J. I. Cirac, and P. ZoUer, Complete Char- 
acterization of a Quantum Process: The Two-Bit Quan- 
tum Gate, Phys. Rev. Lett, vol.78, 390 (1997). [19] 

[14] G. M. D'Ariano and P. Lo Presti, Quantum tomography 
for measuring experimentally the matrix elements of an 
arbitrary quantum operation, Phys. Rev. Lett., vol.86, 
4195 (2001). 

[15] J. B. Ahepeter, D. Branning, E. Jeffrey T. C. Wei, P. G. 

Kwiat, R. T. Thew, J. L. O'Brien, M. A. Nielsen, and A. [20] 
G. White, Ancilla-Assisted Quantum Process Tomogra- 
phy, Phys. Rev. Lett. vol. 90, 193601 (2003). 

[16] G. M. D'Ariano and P. Lo Presti, Imprinting complete 

information about a quantum channel on its output state, [21] 
Phys. Rev. Lett. vol. 91, 047902 (2003). 

[17] M. Mohseni and D. A. Lidar, Direct characterization 
of quantum dynamics, Phys. Rev. Lett, vol.97, 170501 



(2006). 

S. G. Schirmer, A. KoUi, and D. K. L. Oi, Experimen- 
tal Hamiltonian identification for controlled two-level sys- 
tems, Phys. Rev. A 69, 050306(R) (2004) 
S. G. Schirmer and D. K. L. Oi, Two-Qubit Hamiltonian 
Tomography by Baysean Analysis of noisy Data, Phys. 
Rev. A. 80, 022333 (2009); S. G. Schirmer and D. K. L. 
Oi, Quantum System Identification by Bayesian Analysis 
of Noisy Data: Beyond Hamiltonian Tomography, Laser 
Physics 20(5), 1203-1209 (2010) 

Jared H. Cole, Sonia G. Schirmer, Andrew D. Greentree, 
Cameron J. Wellard, Daniel K. Oi, and Lloyd C. HoUen- 
berg, Identifying an experimental two-state Hamiltonian 
to arbitrary accuracy, Phys. Rev. A vol.71, 062312 (2005) 
Jared H. Cole, Andrew D. Greentree, Daniel K. L. Oi, 
Sonia G. Schirmer, Cameron J. Wellard, and Lloyd C. 
L. Hollenberg Identifying a two-state Hamiltonian in the 
presence of decoherence, Phys. Rev. A vol.73, 062333 



11 



(2006) 

[22] M. Zhang, S. G. Schirmcr, H. Y. Dai, W. Zhou and 
M. Lin. Experimental design and idcntifiabihty of model 
parameters for quantum systems [C] Joint Proceeding of 
4.8th IEEE Conference on Decision and Control and 28th 
Chinese Control Conference. 2009: 3827-3832. 

[23] R. Alicki and K. Lendi, Quantum Dynamical Semi- 
groups and Applications, Lecture Notes in Physics No. 
286 (Springer- Verlag, Berlin, 1987). 



[24] W. Feller, An Introduction to Probability Theory and 
Its Applications, Vol.1, John Willey & Sons, 3rd edi- 
tion, 1968 

[25] J.O.Berger, Statistical Decision Theory, Springer- Verlag, 
New York,1980 

[26] A.Wald, Sequential Analysis, John Wiley, New York,1950 
[27] F.J.Anscombe, Sequential Estimation, John Roy. Stat. 
Soc, Series B,15,l-29,1953 



